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We derive the non-linear dynamical equations for an augmented electromagnetic railgun, whose 
augmentation circuit is inductively coupled to the gun circuit. We solve these differential equa- 
tions numerically using example parameter values. We find a complicated interaction between the 
augmentation circuit, gun circuit, and mechanical degrees of freedom, leading to a complicated 
optimization problem. For certain values of parameters, we find that an augmented electromag- 
netic railgun has an armature kinetic energy that is 42% larger than the same railgun with no 
augmentation circuit. Optimizing the parameters may lead to further increase in performance. 



I. INTRODUCTION 

Launch systems, such as the electromagnetic rail- 
gun (EMG), have complex electrical transient phenom- 
ena []], 0|. During launch, the transient involves the 
build-up and penetration of a magnetic field into the 
metallic material. The dynamics of magnetic field pene- 
tration into the metal rails is described by a well-known 
magnetic diffusion equation [3[ , which leads to a skin ef- 
fect where a large current is transported inside a narrow 
channel. In a railgun that has a moving conducting ar- 
mature, the effect is called a velocity skin effect (VSE) 
and is believed to be one of the major problems in limit- 
ing railgun performance [3|, |4| , because it leads to intense 
Joule heating of the conducting materials, such as rails 
and armatures. To what extent the VSE is responsible 
for limiting the performance of solid armatures is still the 
subject of research 

For high-performance EMGs, such as those planne d by 
the navy for nuclear and conventional warships [91-KLlj, 
the armature velocity must be increased while keeping 
the length of the rails fixed. One approach to increas- 
ing the armature velocity is to use some sort of aug- 
mentation to the EMG circuit. Various types of aug- 
mentation circuits have been considered [12|, including 
hard magnet augmentation fields [13| and superconduct- 
ing coils jl3.[l5j. 

In this paper, we consider a lumped circuit model of 
an augmented EMG. The augmented EMG consists of 
two inductively coupled circuits. The first circuit is the 
EMG circuit, containing the rails connected to a voltage 
source, V g (t), that powers the rails and armature. The 
augmentation consists of an external circuit that is in- 
ductively coupled to the EMG circuit via magnetic field, 
but is not connected electrically to the EMG circuit. The 
augmentation circuit has its own voltage generator, V a (t). 
See Figure 1 for a schematic layout. Figure 2 shows the 
equivalent lumped-circuit model that we are considering. 

A simplified arrangement has been previously consid- 
ered, where the EMG circuit was augmented by a con- 
stant external magnetic field [13| . In our work, we assume 
that a real augmentation circuit produces the magnetic 
field that couples to the EMG circuit, and hence, the 
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FIG. 1. Schematic diagram of an inductively augmented EMG 
with its surrounding augmentation circuit. Magnetic flux 
from the augmentation circuit inductively links to the gun 
circuit. 



EMG circuit is interacting with the augmentation circuit 
through mutual inductance, see Figure 1. This results in 
a "back action" on the augmentation circuit by the EMG 
circuit. This back action modifies the current in the aug- 
mentation circuit, resulting in a non-constant 5-field due 
to the augmentation circuit. Furthermore, the EMG cir- 
cuit is coupled to the mechanical degree of freedom, the 
moving armature, which leads to a variation (with ar- 
mature position) of the self inductance and resistance of 
the EMG circuit. The variations of the self inductance, 
resistance and mutual inductance in the EMG circuit act 
back on the augmentation circuit, leading to a complex 
interaction between the gun circuit, the augmentation 
circuit, and the mechanical degree of freedom (the arma- 
ture). The resulting dynamical system is described by 
three non-linear differential equations that are derived 
in Section ITTT1 In order to derive these dynamical equa- 
tions, we first consider Ohm's law in the presence of an 
external magnetic field, which is the field produced by 
the augmentation circuit. 
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II. OHM'S LAW IN A MAGNETIC FIELD 

Consider a series electrical circuit with a voltage gen- 
erator, V, resistor R, and inductor L. Assume that an 
external source of current creates a time-dependent mag- 
netic field, B e , and at some time a quantity of magnetic 
flux from this field, fa, passes through the electrical cir- 
cuit. The time-dependent external magnetic field, B e , 
creates a non-conservative electric field, E e , which are 
related by the Maxwell equation 

curlE e = -^ (1) 

Taking an integration contour along the wire of the elec- 
trical circuit, and integrating Eq. (pQ) over the area en- 
closed by this contour, using Stoke's rule to change the 
surface integral to a line integral, we have the emf in- 
duced in the electrical circuit, s e , given in terms of the 
rate of change of the external flux (due to the external 
magnetic field B e ): 




There are two sources that do work on the electrical cir- 
cuit, the voltage generator, V ' , and the emf, e e , induced in 
the circuit by the external magnetic field, B e . The power 
entering the circuit, e e I + I V, goes into Joule heating in 
the resistor and into increasing the stored magnetic field 
energy: 

£e I + IV = I 2 R+j t hl 2 (3) 

Dividing Eq. (J3j) by /, and solving for the externally in- 
duced emf, e e , in terms of the circuit quantities and the 
rate of change of external flux, fa, 

„ ^ dl d6 e , N 

Equation (|4j) can be rewritten in terms of the total emf 

as 

_y + /i * = _ L ^_^_^ (5 ) 
dt dt dt w 

where —Ldl/dt = dfa/dt is the emf generated in the 
circuit due to the internal current / in the circuit, and 
—d(j) e /dt is the emf generated in the circuit by the ex- 
ternal magnetic field B e . The sum of the internal and 
external flux is the total flux: <fr = fa + fa. Equation (|5]) 
relates the rate of change of total flux through the circuit, 
d(j)/dt, to the voltage drops along the circuit, — V + IR, 
and is an expression of Ohm's law in the presence of a 
magnetic field. We will use this form of Ohm's law to 
write the coupled equations for the augmented railgun. 

III. DYNAMICAL EQUATIONS 

Consider an augmented railgun composed of an aug- 
mentation circuit with voltage generator V a (t) and a gun 
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FIG. 2. The equivalent circuit is shown for the augmented 
railgun in Figure [T] Magnetic flux from the augmentation 
circuit couples to the gun circuit through mutual inductance 
M{x). The self inductance of the gun circuit, L g (x), the mu- 
tual inductance, M(x), and the resistance, R g (x), are func- 
tions of the armature position, x(t). 



circuit with voltage generator, V g (t). We assume that 
the circuits are inductively coupled, but have no electri- 
cal connection, see Figure HJ The equivalent circuit for 
the augmented railgun is shown in Figure [2] The resis- 
tance of the gun circuit, R g {x), changes with armature 
position x(t) and can be written as 

R g (x) = R g0 + R' g x(t) (6) 

where R g o is the resistance of the gun circuit when x = 0, 
and R' is the gradient of resistance of the gun circuit at 
x = 0. 

The total flux in the gun circuit, (j) g , and the total flux 
in the augmentation circuit, fa, can be written as 

fa = Lglg + M g J a (7) 

fa = L a I a + M ag I g (8) 

where I a and I g are the currents in the augmentation 
circuit and gun circuit, respectively, L a and L gi are the 
self inductances of the augmentation and gun circuits, 
respectively, and M ga and M ag are the mutual induc- 
tances, which must be equal, M ga = M ag = M(x). The 
self inductance of the gun circuit, L g (x), changes with ar- 
mature position x(t). Also, the area enclosed by the gun 
circuit changes with armature position, and therefore, the 
coupling between the augmented circuit and gun circuit, 
represented by the mutual inductance, M(x), changes 
with armature position, x(t). Furthermore, in order for 
the free energy of the system to be positive, the self in- 
ductances and the mutual inductance must satisfy [HI 

M{x) = k^L a L g (x) (9) 

for all values of x. Here, the coupling coefficient must 
satisfy |fc| < 1. We can write the self inductance of the 
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gun circuit as 

L g (x)=L g0 + L g x(t) (10) 

where L g o is the inductance when x = and L' g is the 
inductance gradient of the gun circuit. Similarly, the 
mutual inductance between augmented circuit and gun 
circuit can be written as 



M(x) = M + M' x(t) 



(11) 



where Mq is the mutual inductance when x = and M' 
is the mutual inductance gradient. For x = and x = £, 
where £ is the rail length, Eq. (|9j) gives 



M k^LaLgo 



fcjLa (L K0 + L' g e)-M 



(12) 
(13) 



The coupling coefficient, fc, can be positive or nega- 
tive, and as mentioned above, must satisfy |fc| < 1. The 



sign of k determines the phase of the inductive coupling 
between augmentation and gun circuits. Choosing the 
coupling coefficient fc, and the two self inductances, L g o 
and L a , Eq. ([T2]) then determines the value of the mu- 
tual inductance, Mo. Then, choosing a value for the rail 
length i, and the self inductance gradient, U , Eq. (p~3|) 
determines the mutual inductance gradient, M' . See Ta- 
ble [J for parameter values used. 

Two dynamical equations for the augmented railgun 
are obtained by applying Ohm's law in Eq. ([5]) to the 
gun circuit and to the augmented circuit: 



-V g + I g R g (x) 



j t [L a I a + (M + M f x(t))I g ] 



L' g X(t))Ig 

(M + M'x(t))I a ] 



(14) 



The third dynamical equation is obtained from the 
coupling of the electrical and mechanical degrees of 
freedom[2]. Therefore, the three non-linear coupled dy- 
namical equations for I g (t), I a (t) and x(t) are given by: 



V a (t) + I a (t)Ra = ~L a ^ - (M + M'x(t)) ^ - M'I g (t)^ 



-Vg(t)+Ig(t) (RgO + R' g X(t)) 

d 2 x(t) 



din 



dx(t) 



(15) 



■L> BX ®)-±-L> g IS) dt 



ni- 



di 2 



(M + M'x(t))^-M'I a (t)^ (16) 

(17) 



where V g (t) and V a (t) are the voltage generators that 
drive the gun and augmentation circuits. From Eq. (fT7|) . 
we see that the EMG armature has a positive accelera- 
tion independent of whether the gun voltage generator 
is a.c. or d.c. because the armature acceleration is pro- 
portional to I g 2 (t). The armature velocity is essentially 
the integral of I g 2 (t), and therefore a higher final velocity 
will be achieved for d.c. current, and associated d.c. gun 
voltage, V g (t) = Vgo, where V g o is a constant. Of course, 
the actual current will not be constant in the gun circuit 
because of the coupling to the moving armature and to 
the augmentation circuit. We want to search for solu- 
tions where the armature velocity is higher for a EMG 
with the augmented circuit than for an EMG without an 
augmentation circuit. In order to increase the coupling 
between the augmentation circuit and the gun circuit, 
we choose an a.c. voltage generator in the augmentation 
circuit: 

V a (t) = KosinfW) (18) 

where V a o is a constant amplitude for the voltage gener- 
ator and uo is the angular frequency of the augmentation 
circuit voltage generator. 

As an example of the complicated coupling between 
augmentation circuit and gun circuits, we will also ob- 



tain solutions for a d.c. voltage generator for the aug- 
mentation circuit. As we will see, the gun circuit causes 
a back action on the augmentation circuit, leading to a 
non-constant current in the augmentation circuit. 

We need to choose initial conditions at time t = 0. We 
assume that there is no initial current in the gun and 
augmentation circuits and that the initial position and 
velocity of the armature are zero: 

Ig(0) = (19) 

I a (0) = (20) 

x(0) = (21) 

dx(0) , , 

-li=0 (22) 

For the special case when L' g = 0, M' = 0, and R g = 0, 
the mechanical degree of freedom described by Eq. (fT7|) 
decouples from Eqs. (fT5]) - (fT7|) . In this case, Eqs ([T5]) - 
(fTB|) describe a transformer with primary and secondary 
circuits having voltage generators, V a (t) and V g (i), re- 
spectively. The solution for the mechanical degree of 
freedom is then x(t) = for all time t. 

In what follows, we solve the dynamical Eqs. ([T5]) -([T7 |) 
numerically for the case when the EMG and augmented 
circuits are coupled. 
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TABLE I. Electromagnetic gun and augmentation circuit parameters. 



Quantity 


Symbol 


Value 


length of rails (gun length) 


t 


10.0 m 


mass of armature 


in 


20 kg 


coupling coefficient 


k 


±0.80 


self inductance of rails at x — 


LgO 


6.0xl0" 5 H 


self inductance of augmentation circuit 


La 


6.0xl0" 5 H 


self inductance gradient of rails 


L' B 


0.60xl0" 6 H/ 


resistance of augmentation circuit 


Ra 


0.10 Q 


resistance of gun circuit at x = 


RgO 


0.10 Q 


resistance gradient of gun circuit at x = 


Rg 


0.002 tt/m 


voltage generator amplitude in gun circuit 


V g0 


8.0 xlO 5 Volt 


voltage generator amplitude in augmentation circuit 


V a0 


8.0 xlO 5 Volt 




t [10" 3 sec] t [10~ 3 sec] 

FIG. 3. For a.c. augmentation circuit voltage, with uj — 10 3 rad/sec, the armature position, velocity, and augmentation and 
gun circuit currents are plotted for two different values of coupling constant, k — +0.80 and k — —0.80. 



IV. NUMERICAL SOLUTION 

As described above, if we choose the parameters /c, L a , 
L g o, then Mo is determined from Eq. (fl2j) . Next, if we 
choose L' g and i, then the inductance gradient, M', is 
determined by Eq. ([T3]) . See Table U for values of param- 
eters used in the calculations below. 

The sign of the coupling coefficient, /c, affects the in- 
teraction of the augmentation and gun circuits in subtle 
ways. We solve Eqs. ([T5]) - ([TT)) for a positive and negative 
value of the coupling coefficient k. Figure [3] shows a plot 
of armature position, x(t), velocity v(t), and currents 



I g (t) and I a (i), for a gun circuit powered by a d.c. volt- 
age V g o, and an augmentation circuit with an a.c. voltage 
source given by Eq. ([T8]) , using the parameters in Table [I] 
and coupling coefficient values k = +0.80 and k = —0.80. 
Note that the rail length is taken to be t = 10 m, so the 
plots are only valid for time < t < tf , where tf is given 
by x(tf) = £. For k = +0.80 we have tf = 5.563 ms, 
and for k = -0.80 we have tf = 5.277 ms, see Table HI 
From Figure [3j for k = +0.80, for short times t ~ 2 
ms the armature velocity is smaller than for k = —0.80. 
However, for longer times t ~ 5 ms, the armature veloc- 
ity is higher for k = +0.80. For k = +0.80, for short 
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FIG. 4. For a.c. augmentation circuit voltage, with uj = 10 3 rad/sec, the left plot shows the armature velocity vs. time for 
three values of coupling constant, k = —0.80, k = +0.80, and k = 0.0, for parameter values in Table [T] Note that the velocity 
curves in this figure for k = —0.80 and k = +0.80, are the same as the curves in Figure [3] The right side plot shows the current 
in the gun circuit and in the decoupled (k = 0) augmentation circuit. The corresponding kinetic energy plots are shown in 
Figure \S\ 




1 2 3 4 5 6 



t [lCT 3 sec] 

FIG. 5. For a.c. augmentation circuit voltage, with uj = 10 3 
rad/sec, the armature kinetic energy is plotted vs. time for 
three values of coupling coefficient, k = —0.80, 0.0, and 0.80, 
for parameters values shown in Table |U The kinetic energy is 
largest when k — +0.80. For the corresponding velocity plots, 
see Figure HI 



time t ~ 2 ms, the augmentation circuit current I a (i) 
is out of phase with the gun current I g (t). However, for 
k = —0.80, for short time t « 2 ms, the augmentation cir- 
cuit current and gun circuit currents are approximately 
in phase. Even though the gun circuit has a d.c. volt- 
age generator, the gun circuit current varies due to the 
coupling with the augmentation circuit and with the me- 
chanical degree of freedom x(t). These plots illustrate the 
complicated coupling between the augmentation circuit, 
the gun circuit, and the mechanical degree of freedom, 
x(t). When the armature reaches the end of the rails at 
i = 10 m, its velocity is v(t f ) = 4769 m/s for k = +0.80, 
whereas its velocity is v(tf) = 3336 m/s for k = —0.80. 
The armature kinetic energy for k = 0.0 is 159.6 MJ, 



TABLE II. A.C. augmentation. For parameters in Table [I] 
taking the angular frequency of the augmentation circuit volt- 
age in Eq. ([18[) as uj — 10 3 rad/sec, the shot time, £/, armature 
velocity, v(tf), and armature kinetic energy (KE) in units of 
mega Joule, are shown for three different values of coupling 
constant k. 



k 


t f [10- 3 s] 


v(t f ) [10 3 m/s] 


KE [M J] 


-0.80 


5.277 


3.336 


111.3 


0.0 


5.545 


3.994 


159.6 


+0.80 


5.563 


4.769 


227.4 



TABLE III. D.C. augmentation. For parameters in Table [H 
taking the angular frequency of the augmentation circuit volt- 
age in Ea. (|18[) as uj = 0, the shot time, £/, armature velocity, 
v(tf), and armature kinetic energy (KE) in units of mega 
Joule, are shown for three different values of coupling con- 
stant k. 



k 


t f [10- 3 s] 


v(t f ) [10 3 m/s] 


KE [M J] 


-0.80 


4.860 


4.056 


164.5 


0.0 


5.545 


3.994 


159.6 


+0.80 


6.182 


3.868 


149.6 



whereas for k = +0.80 the armature kinetc energy is 
227.4 MJ, so there is a 42% increase in kinetic energy 
in the augmented EMG compared to the non-agumented 
EMG. 

We can decouple the augmentation circuit from the 
gun circuit by setting k = 0, and through Eq. ([T2]) and 
([T3]h this leads to zero values for mutual the inductances, 
M = and M' 0, in Eq. (fT5|) (fT7]). In this case, 
Eq. (US) is decoupled from Eq. (J6]) and Eq. ([17]). Equa- 
tion JTSJ) then describes a simple L-R circuit with a volt- 
age generator V a (t), leading to a current in the (decou- 
pled) augmentation circuit that is purely sinusoidal. In 
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01234567 01234567 
t [10~ 3 sec] t [10 -3 sec] 



FIG. 6. For d.c. augmentation circuit voltage (with uj = 0), the armature position, velocity, and currents are shown for the 
case when the voltage generators in augmentation circuit, V a (t), and gun circuits, V g (t), are d.c. generators (cj = 0) having 
values shown in Table |U 



this case, Equations (fT6j) and (fT7|) describe a simple EMG 
that is not coupled to an augmentation circuit. Figure |4] 
shows the velocity vs. time for this EMG with no aug- 
mentation circuit. Also shown are the currents in the 
gun circuit and (decoupled) augmentation circuit. When 
k = 0, the time for the armature to reach the end of the 
rails is tf = 5.545 ms and its velocity is v = 3994 m/s. 
For the parameters in Table HJ the augmented EMG with 
k = 0.80 has a 19 % increase in velocity over the k = 
non-augmented EMG. As mentioned above, in terms of 
kinetic energy, the augmented EMG with k = 0.80 has a 
42% larger kinetic energy than the non-augmented gun 
(with k = 0.0), see Figure [5j where we plot the kinetic 
energy of the armature as a function of time. 

Finally, we show the effect of the coupling between aug- 
mented and gun circuits when the augmented circuit has 
a d.c. voltage generator with constant amplitude volt- 
age, V a (t) = V a0 . We solve Eqs. (jl5j)-([I7j) for coupling 
coefficient values k = +0.80 and k = —0.80. For d.c. 
augmentation (with uj = 0), for k = +0.80, the arma- 
ture velocity at tf = 6.182 ms is v = 3868 m/s, while 
for k = —0.80 the armature velocity at tf = 4.860 ms is 
v = 4056 m/s, see Table HIB Therefore, for a d.c. volt- 
age generator in the augmentation circuit, the negative 
value k = —0.80 leads to a higher armature velocity than 



the positive k value, which is just opposite for the case 
of an a.c. voltage generator in the augmentation circuit. 
The plots are shown in Figure [6] and the kinetic energy 
is shown in Figure [71 Note that even though the volt- 
age generators in the gun and augmentation circuits are 
now taken to be d.c, the resulting currents, I g (t) and 
I a (i), are not constant, due to the interaction between 
the augmentation and gun circuits and these circuits in- 
teracting with the mechanical degree of freedom x(t) of 
the armature. 



V. SUMMARY 

We have written down the non-linear coupled differen- 
tial equations for an augmented EMG that is inductively 
coupled to the augmentation circuit. We solved these 
differential equations numerically using example values 
of parameters. For the sample parameters that we used, 
with k = +0.80 the a.c. augmented EMG had a 42% 
larger kinetic energy of the armature than the same EMG 
with no augmentation. We have made no effort to vary 
the parameter values to optimize the augmented EMG 
performance. The improvement in performance may be 
higher for other parameter values. In this work we have 
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neglected many practical design concerns, such as heating 
and melting of the rails [2] . We have considered the sim- 
plest configuration of an inductively coupled augmented 
EMG, with a single external circuit coupling to the EMG 
circuit. Other configurations may be explored, such as 
multi-stage inductively coupled circuits to the EMG cir- 
cuit. Further work is needed to explore in detail the 
parameter space of this augmentation scheme as well as 
alternative augmentation schemes. 



FIG. 7. For d.c. augmentation circuit voltage (with uj = 0), 
the armature kinetic energy is plotted vs. time for three values 
of coupling coefficient, k = —0.80, 0.0, and 0.80, for param- 
eters values shown in Table U The kinetic energy is largest 
when k = —0.80. For the corresponding velocity plots, see 
Figure [6] 
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